Charge response function and a novel plasmon mode in graphene 
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Polarizability of non-interacting 2D Dirac electrons has a Xj y/qv — to singularity at the bound- 
ary of electron-hole excitations. The screening of this singularity by long-range electron-electron 
interactions is usually treated within the random phase approximation. The latter is exact only 
in the limit of A*' — > oo, where A'^ is the "color" degeneracy. We find that the ladder- type vertex 
corrections become crucial close to the threshold as the ratio of the n-th order ladder term to the 
same order RPA contribution is In" \qv — uj\/N^ . We perform analytical summation of the infinite 
series of ladder diagrams which describe excitonic effect. Beyond the threshold, qv > uj, the real 
part of the polarization operator is found to be positive leading to the appearance of a strong and 
narrow plasmon resonance. 

PACS numbers: 73.23.-b, 72.30.-|-q 
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Introduction. Many properties of interacting two- 
dimensional electrons with linear Dirac spectrum e = 
±wp, found in graphene monolayers, differ sharply from 
those with the parabolic spectrum present in conven- 
tional semiconductor heterostructures [l|, Q ■ Vanishing 
density of states at the Dirac point and less effective 
screening of Coulomb interaction lead to electronic prop- 
erties of graphene being qualitatively different. 

One of the distinct signatures of non-interacting 2D 
electrons in graphene, caused by the absence of spec- 
trum curvature, is a divergent behavior of the polar- 
ization operator (charge susceptibility) at the threshold 
for the excitation of electron- hole pairs [1, 0j S, B 01' 
n(")(tj,(7) = -Nq^ /l&^q^v'^ - u'^ , where N is the de- 
gree of degeneracy (in graphene iV = 4 in the absence 
of magnetic field). The effects of electron-electron inter- 
action on the polarization operator are customarily ac- 
counted for by the random phase approximation (RPA) 
[Ej B B which sums the infinite series of electron 
loops, 
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where Vq = 27re^/(j, dielectric constant of the substrate 
being incorporated into the charge e^. In particular, RPA 
predicts that the imaginary part of the polarization op- 
erator (which determines absorption in the system and 
the density-density correlation function) instead of di- 
verging actually vanishes at the threshold according to, 
— n^p^ oc g~'^^uj — qv, where g = e^/v is the dimen- 
sionless interaction constant. 

A notable property of 2D electrons, well captured by 
RPA in a conventional Fermi liquid, is the existence 
of a low-frequency collective mode of charge oscillations 
whose spectrum is, = 2e^Epq, where Ep is the Fermi 
energy. In undoped graphene, Ep = 0, the plasmon 
mode is thus absent within RPA (finite doping or tem- 
perature can lead to the usual RPA plasmons 11, 1^). 
Mathematically this comes from the fact that the real 
part of n^"^ is negative for qv > uj (i.e. beyond the do- 
main of electron- hole excitations). In the present paper 



we demonstrate that this conclusion is an artifact of the 
approximation ([T]) and that a collective mode does exist 
in undoped graphene. This mode describes charge fluc- 
tuations in a system of electron-hole pairs interacting in 
a final state ("excitonic" effect), rather than free pairs, 
as implied by Eq. 

We begin with noting that RPA is much less justified 
in case of undoped graphene than in the case of conven- 
tional parabolic spectrum, where it is formally valid for 
momenta less than the inverse screening length, g <C /c. 
In undoped graphene the vanishing density of carriers 
ensures that k ^ 0. In particular, direct calculation of 
the electron self-energy to the second order in interaction 
demonstrates that non-RPA contributions are generally 
smaller than the RPA terms only by virtue of 1 /N [iSj . 
The RPA is exact only in the formal limit of ^ cxd, 
since each RPA term corresponds to the largest number 
of loops possible in each order of the perturbation theory. 

Below we analyze consecutive orders of the pertur- 
bation series in the bare interaction Vq = 2ne'^ /q and 
compare RPA diagrams to the ladder corrections. We 
find that the ratio of the n-th ladder term to the 
corresponding RPA contribution of the same order is 
iV-"ln"(T-2^) and large close to the threshold. Af- 
ter extracting the leading singular terms in each or- 
der, we calculate the infinite ladder series and obtain a 
non-perturbative expression for the polarization operator 
near the threshold. Summation of the leading divergen- 
cies is possible once we realize that for qv / \qv — u}\ 3> 1 
the main contribution arises from the processes that in- 
volve quasi-one-dimensional motion of interacting parti- 
cles along the direction of the external momentum q. In 



general, when A^ ~ In ( 



r), the perturbation expan- 



sion of the polarization operator receives similar contri- 
butions from both the ladder and RPA terms. Thus, the 
polarization operator is composed of all possible combi- 
nations of ladders and loops. 

Non-interacting electrons. To demonstrate our method 
and to introduce some notations it is instructive to begin 
with reproducing the zeroth-order polarization operator. 
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Close to the resonance, when \uj — qv\ <C qv, only those 
transitions are important that have momenta p and q 
directed almost along the same line. More precisely, the 
electron is taken from the state in the lower cone with 
the momentum p almost antiparallel to q and placed in 
the upper cone in the state with momentum p + q which 
is almost parallel to q. Obviously, this is possible only 
as long as p < q. The relevance of these particular tran- 
sitions is understood from the form of the polarization 
operator. 



n(°)(w,g) = -iNTi 
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where the electron Green's function consists of the con- 
tributions from both cones /3 = ±1, 



G, 
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where ffp = <t • p/p is the projection of the pseudospin 
Pauli matrix onto the direction of electron momentum. 
The energy integral in Eq. ^ yields, 

nWr. „^-^W Tr(l + /?gp+q)(l-/3gp) 

4^ y ^ ^ - I3v{p + \P + ci\) + ^|3v' ^' 

we observe that i) for positive uj the singular denom- 
inator appears for /? = -fl, ii) the trace operation in 
the numerator imposes that p and p -I- q are antipar- 
allel (for the parallel configuration the trace vanishes). 
Close to the threshold (w — qv) collinear processes be- 
come dominant allowing us to approximate |p -I- q| ~ 
Q ~P + 2fq-p) ' 'where 6 denotes the angle between p and 
— q and Tr(l -|- (Tp+q)(l — (jp) « 4, thus we obtain from 
Eq. Q, 
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Upon calculating the integrals we recover II'-'^) (w, g) 



N 



167211 1^^"^ I ~ ^- Note that, as expected, the rele- 
vant angles are small 9 ^ ^|1 — ijj/qv\ ^ 1. 

First-order correction. To the first order in the 
electron-electron interaction there are three diagrams. 
Fig. 1. Of these diagrams, the self-energy correction. 
Fig. la), is the most singular near the threshold. 



Ii%{u,q) = - 



^5/2 



162\/2 {qv - w)3/2 



ln(/C/g), (6) 



where /C is the upper momentum cut-off. This singularity 
indicates that the infinite summation of the self-energy 
diagrams have to be performed first, thus yielding the 
renormalization of the electron spectrum. To the low- 
est order in interaction this gives Vp = v[l + jln {)C/p)] 




FIG. 1: First-order interaction correction: self-energy correc- 
tion a), sub-leading RPA correction b), leading vertex cor- 
rection c). Figure d) illustrates the origin of the singular 
behavior, ~ {qv — uj)~^ \n\qv — from almost collinear elec- 
tron propagation. Both states pi, p2 belong to the lower (— ) 
cone, while pi + q, P2 -I- q correspond to the upper (-I-) cone. 



In the case of a conventional Fermi liquid the cor- 
responding velocity renormalization is trivial. In case 
of a linear spectrum considered here, however, logarith- 
mic terms lead to a slight curvature of electron spectrum 
which can smear the threshold singularities. In what fol- 
lows we assume the renormalized velocity but neglect this 
smearing by taking the velocity at a typical momentum, 
determined by the external momentum v = Vq. 

The next term, lb), is the RPA-corrcction which can 
be readily written as. 



n 



(1) 

RPA 



nN'^g q^ 
16^ qv — 



- LO 



(7) 



It can now be verified that the vertex correction. 
Fig. Ic), yields a contribution which is opposite in sign 
and more singular than Eq. ([7]) as a result of the long- 
range character of Coulomb interaction. The energy in- 
tegrations are performed independently in each half of 
the diagram, resulting in 
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To extract the leading contribution wc note that it comes 
from almost collinear processes with a = /3 = +1 and 
both p -I- q, p' + q are almost parallel to q and almost an- 
tiparallel to both p, p'. Using the small-angle expansion 
introduced in the previous section, and approximating 
Vp-p' ~ /\p — p'\ we obtain 
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The integrals over the two angles are independent and 
yield the same resonant denominator as in Eq. ([7]) . How- 
ever, the integrals over the magnitude of electron mo- 
menta contain a logarithmic divergence which is due to 
the long-range nature of Coulomb interaction. To cut- 
off this divergence we note that, |pi — P2I ~ bi — P2I + 
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i^'"^^ "i , where 9i2 is the angle between two momenta. 

2|pi-P2| ' -^^ ° 

RecaUing that relevant angles ~ l^w/w — 1| and that 
typical Pi ~ q, we observe that small-angle expansion 
fails for \pi — P2I ~ g-y/jgu/w — 1|, which should be used 
as the lower cut-off. As a result we obtain, 
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We, therefore, find that RPA loops exceed the vertex 
terms only when N ^ In ( ^^'^^^^ )■ In the opposite limit, 
close to the resonance, N <C In ( , '^^ , ) the ladder dia- 

' ^ \uj — qv\ ' 

grams dominate. We now turn to the summation of the 
ladders series. 

Summation of the infinite ladder. In the second 
order in electron-electron interaction the most domi- 
nant contribution again comes from the ladder diagram. 
Fig. 2a), with three pairs of electron lines each yield- 
ing 1/^/qv — w-singularity, and two Coulomb interactions 
providing two additional powers of the logarithm. The 
other diagrams are less singular at to = qv. Diagram 
2b) requires that at least one of the interaction lines car- 
ries large momentum, ~ q, while both momenta in the 
ladder-type diagram 2a) are small. Exchange-type con- 
tributions, 2c)-e), though providing the main power-law 
singularity, involve transferred momenta which are ~ q 
and lack the additional logarithms (while also lacking ex- 
tra factors A^). Repeating the arguments from the pre- 



a) b) 

<TT> <ZX:^XZ> 



c) 



d) 



e) 



P2+Q 



P1 



Pn.1 



FIG. 2: Leading, a), and subleading, b)-e), non-RPA dia- 
grams of the second order. The ladder diagram f) explains 
the origin of the most singular, {qv — uj)~" In" \qv — u>\ contri- 
bution of the arbitrary order n. All states propagating along 
the upper/lower part of the ladder belong to the upper/lower 
cone. 

ceding section we realize that to each order in interaction 
the leading contribution comes from the ladder diagram. 
Fig. 2f), 



„ ji+i , 

4")(c.,q) = -A^Try [nl^rfp. 



^l-P2 ■■■^n-Pn+l 



All the momenta are almost collinear, with those propa- 
gating along the lower part of the ladder, pi, being an- 
tiparallel to external wavevector q and to all the mo- 
menta Pi -l- q along the upper part of the ladder. Energy 
integrations are independent in each step and result in 
the singularity to be the strongest when all the propa- 
gators in upper/lower parts of the ladder belong to the 
upper (-|-)/lower(-) cones. The subsequent trace opera- 
tion is easily performed as 

Tr (1 - CTpJ...(l - Crp„ + J (1 + CTp„+i+q)... 

...X(l + ap,+q)«4"+l, 

yielding the n-th order contribution to the polarization 
operator in the form 

4"^(^,q) = {-irN I v^p,_p,...v_p„^, 

dpi 
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Integrating over the angles and then evaluating the inte- 
grals over the absolute values of momenta in the leading 
logarithmic approximation we obtain. 
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^2qv{qv -cj)T (^) 
where we introduced the following dimensionless variable 
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nite series can now be readily performed. 



The summation of the infi- 



n(c^,<z) = -Nq 
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The above result p3|) for the polarization ladder is 
sufficient to calculate the polarization operator in a more 
general case when N ~ ln( ) ^1. In this regime 
both vertex and loop diagrams have to be summed simul- 
taneously. Such summation can be performed by realiz- 
ing that it yields a geometric series similar to the usual 
RPA one ^ but with the bare polarization operator II'^*'^ 
replaced by the polarization ladder Uy given by Eq. ([13]), 
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Collective mode. Non-zero imaginary part of the po- 
larization operator determines those values of w and q for 
which the dissipation of external field is possible. Imagi- 
nary part of Eq. (flSjl . and as a result of Eq. ^4]) . arises 
in two cases: i) Imaginary x corresponds to the conven- 
tional domain, to > qv, but modified by interaction of 
electrons and holes, ii) Positive 1 < x < 00 correspond to 
the absorption below the threshold and represent a non- 
perturbative effect which arises from lowering the energy 
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FIG. 3; (Color online) Imaginary part of the polarization op- 
erator in units of q/A for = 4 and g — 0.3: non-interacting 
value Ho (black dashed line), RPA value Il'j^p^ from (blue 
dot-dashed line) and V" from Eq. (|14p that describes collec- 
tive response of excitons (red solid line) . The inset illustrates 
the relative height of the excitonic plasmon line to the back- 
ground of individual electron-hole absorption. The position of 
plasmon qv = l.OOGw, the absorption threshold qv — l.OSloj. 



Conclusion. The usual RPA formalism adequately de- 
scribes dynamical charge fluctuations only in the limit 
In I I . Close to the Dirac cone we find that 

\qv-uj\ 

the n-th order term in the ladder series is larger than 
the corresponding term in the RPA series by a factor 
of lii^ {qv/\qv — uj\)/N". This stronger singularity in the 
ladder series arises due to the long range interaction be- 
tween electrons which move almost collinearly to the ex- 
ternal momentum q. The series is summed up analyti- 
cally yielding a non-perturbative result: the density and 
spin response functions acquire non-zero imaginary part 
in the additional frequency range, qu < u < qv. This 
extension is a manifestation of the excitonic effect. The 
reversal of the sign of the electron polarizability in this 
new domain gives rise to a sharp plasmonic mode which 
is absent in the conventional RPA. Interestingly, a similar 
sign reversal of the polarizability due to interactions has 
been reported to yield a surface plasmon in a different 
case of a two-dimensional Anderson insulator [151 1. 



of electron-hole excitations with finite q from attractive 
Coulomb interaction. The new threshold of absorption 
uj = qu is determined from the condition x = I which for 
weak interactions g ^ I, gives 



u = «[l-(glng)V2]. 



(15) 



Even more striking property of Eq. (|13p is the behavior 
of its real part which describes the polarizability of an 
interacting electron-hole pair. It is easy to see that the 
real part is positive for a; > 1 



14|. This change in sign, 



as compared with the real part of the RPA polarizability 
([1]), can already be traced to the first order vertex cor- 
rection (fTO|) , which is opposite to the first order RPA cor- 
rection ^ . This fact has immediate consequence for the 
collective response of excitons, given by Eq. ([13]), which 
develops a pole for a certain value of frequency uj{q). This 
new mode describes the propagation of coherent oscillat- 
ing charge density in a system of interacting electron-hole 
pairs. The corresponding peak in the absorption has a 
finite but small width which is due to the fact that the 
collective mode falls within the range of decay via indi- 
vidual electron- hole pairs. Its behavior is illustrated in 
Fig. 3 for A 4. 

Interestingly, the width of the plasmon is maximal 
around g ~ 0.3. For smaller values of the coupling 
strength the peak narrows by virtue of the narrowing 
of the whole exciton domain, Eq. (|15p . For larger val- 
ues, g 1, the peak width decreases as oc while 
the height increases only as oc g. The position of the 
peak in this limit tends to the value determined from 
the equation (recovered from (|13m4p when x oo), 
= which yields. 



Useful discussions with A. Chubukov, M. Raikh and 
O. Starykh are gratefully acknowledged. This work was 
supported by DOE, Grant No. DE-FG02-06ER46313. 



[1] 
[2] 

[3] 
[4] 

[5] 

[6] 

[7 

[9: 
[lo: 

[11 

[12 
[13 
[14 



Nucl. 



Phys. 



B 73, 



[15] 



uj{q) = qv{l 



(16) 



M. Wilson, Phys. Today 59, 21 (2006). 
A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, A. K. Geim, arXiv:0709.1163 
K. W. Shung, Phys. Rev. B 34, 1264 (1986). 
J. Gonzalez, F. Guinea, and M.A.H. Vozmediano, 
Phys. B 424, 595 (1994). 

J. Gonzalez, F. Guinea, and M.A.H. Vozmediano, 
Rev. B 59, R2474 (1999). 

A.V. Chubukov and A.M. Tsvelik, Phys. Rev. 
220503(R) (2006). 

D.V. Khveshchenko, Phys. Rev. B 74, 161402(R) (2006). 

D. V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001). 
S. Das Sarma, E.H. Hwang, and W-K. Tse, Phys. Rev. 
B 75, 121406(R) (2007). 

Y. Barlas, T. Pereg-Barnea, M. Polini, R. Asgari, and A. 
H. MacDonald, Phys. Rev. Lett. 98, 236601 (2007). 

E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 
(2007). 

O. Vafek, Phys. Rev. Lett. 97, 266406 (2006). 
E.G. Mishchenko, Phys. Rev. Lett. 98, 216801 (2007). 
For a:: > 1 the analytical continuation should be per- 
formed from X G [0,1) to (l,oo). From the require- 
ment that n" is negative we observe that it should yield 
\/l — — ^ —i\/x'^ ~ 1. This implies that the analytical 
continuation must be performed in the upper half-plane 
of complex x, so that arcsin x — > 7r/2-|-i In {x + \J x'^ — 1). 
T.V. Shahbazyan and M.E. Raikh, Phys. Rev. B 53, 7299 
(1996). 



